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ABSTRACT 

Objectives: To assess whether residential proximity 
to industrial incinerators in England is associated with 
increased risk of cancer incidence and mortality. 
Design: Retrospective study using matched case- 
control areas. 

Setting: Five circular regions of radius 10 km near 
industrial incinerators in England (case regions) and 
five matched control regions, 1998-2008. 
Participants: All cases of diseases of interest within 
the circular areas. 

Primary and secondary outcome measures: 

Counts of childhood cancer incidence (<1 5 years); 
childhood leukaemia incidence (<15 years); leukaemia 
incidence; liver cancer incidence; lung cancer 
incidence; non-Hodgkin's lymphoma incidence; all- 
cause mortality; infant mortality (<1 year) and liver 
cancer mortality. 

Results: The estimated relative risks for case circles 
versus control circles for the nine outcomes 
considered range from 0.94 to 1.14, and show neither 
elevated risk in case circles compared to control areas 
nor elevated risk with proximity to incinerators within 
case circles. 

Conclusions: This study applies statistical methods 
for analysing spatially referenced health outcome data 
in regions with a hypothesised exposure relative to 
matched regions with no such exposure. There is no 
evidence of elevated risk of cancer incidence or 
mortality in the vicinity of large industrial incinerators 
in England. 
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INTRODUCTION 

The disposal of solid waste is an increasing 
problem in developed countries as space for 
landfill sites runs out. Incineration is an 
attractive solution as it reduces the weight 
and volume of waste, and prevents growth of 
pathogenic bacteria and the spread of 
disease vectors, while also providing an 
opportunity to recover thermal energy for 
electricity generation/ ^ 



Article focus 

■ Research specific to waste incinerators has pro- 
vided mixed evidence for the effects of proximity 
to incinerators on health. 

■ Older incinerators have been associated with 
increased incidence and mortality from selected 
cancers, while more recent reports show little 
association. 

■ Despite this, the effect of incinerator emissions 
remains a public health issue, 

Key messages 

■ Our results show no evidence of elevated risk for 
individuals living in areas containing an inciner- 
ator compared to individuals living in matched 
areas without an incinerator. 

■ Moreover, within areas, there is little evidence of 
an increase in risk for those living in close prox- 
imity to an incinerator compared to those living 
further away. 

Strengths and limitations of this study i 

■ The strengths of this study include adjustmen 
for the age/sex structure and deprivation of the 
underlying population at a small-area level. The 
design using matched areas reduces the need to 
adjust for circle-level risk factors for which data 
are unavailable. 

■ Limitations include the long lead times for adult 
cancers and use of distance rather than modelled 
or directly measured pollution levels as a proxy 
for exposure. Data are at small-area level, making 
adjustment for individual risk factors such as 
smoking status infeasible. 



Incineration produces emissions of carbon 
dioxide (CO2) and water (H2O). Emissions 
can also include acid gases (hydrogen chlor- 
ide, hydrogen fluoride and sulfur dioxide), 
nitrogen oxides, heavy metals (cadmium, 
mercury, arsenic, vanadium, chromium, 
cobalt, copper, lead, manganese, nickel, thal- 
lium and tin), polychlorinated biphenyls 
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(PCBs), polycyclic aromatic hydrocarbons (PAHs) and 
dioxins, although emissions of dioxins can be minimised 
by keeping the combustion temperature above 1100°C. 
Residual ash from solids is also produced, and may 
contain aluminium, calcium, silicates, iron, sodium, mag- 
nesium, potassium, lead and zinc/ The International 
Agency for Research in Cancer (lARC) considers some 
of these substances to be potentially carcinogenic to 
humans^; in particular, some heavy metals and PAHs are 
classified as either Group 1 (carcinogenic to humans) or 
Group 2A (probably carcinogenic to humans) with 
several others classified as Group 2B (possibly carcino- 
genic to humans). PCBs are classified as Group 2A. 
Exposure to these substances may occur either by air or 
by consumption of contaminated food. Possible effects 
include reproductive impairment and carcinogenesis.^ 

There is ongoing public concern over the extent to 
which emissions produced by industrial incinerators may 
have an adverse effect on human health. Several studies 
have investigated this potential risk by applying a variety 
of spatial statistical methods to both small area-level and 
individual-level data,^~^ although they have had limited 
ability to control for confounding factors such as depriv- 
ation and other industrial emissions. 

While studies of older incinerators have found evi- 
dence of an adverse health effect, reports based on data 
from newer incinerators show little association. Recent 
reviews have noted that there is low risk of adverse 
health effects due to residency close to incinerators and 
any effects are probably not measurable.^ ^ They base 
this view on assessments of the effects of air pollutants 
on health and the fact that modern (post-European 
Directive 2000/76/EC) incinerators only contribute a 
small amount to local air pollution. 

Though there is now considerable literature investigat- 
ing the possible elevation in disease risk in the vicinity of 
an industrial point source, work specific to incinerators 
has provided mixed evidence for the effects of proximity 
to incinerators on health. Elliott et af found excess 
cancers in people living in close proximity to municipal 
solid waste incinerators (observed/ expected ratios for 0- 
3 km from the incinerator ranged from 1.04 for all cancers 
and colorectal cancer to 1.13 for liver cancer), but sus- 
pected residual confounding. Kim et af^ also found 
adverse health effects in the vicinity of solid waste incinera- 
tors, with an attributable burden of 0.2% and 0.1% of the 
total burden of respiratory disease and cardiovascular 
disease, respectively, in Seoul in 2007. Porta et at^ report 
an increased risk in cancer of 3.5% for people living close 
to old incinerators but note that studies suffer from limita- 
tions due to poor exposure assessment, ecological level of 
analysis and lack of information on relevant confounders. 
Most of these, and other studies, ^^^^ are concerned with 
older incinerators, and their results may not be applicable 
to incinerators currently in operation and under tighter 
regulatory control. 

Other studies found no evidence of association of 
health risks with incinerators: Gouveia and Ruscitto^ 



investigated various cancers in Brazil; Vilavert et al 
monitored the concentrations of polychlorinated 
dibenzo-/^dioxins and dibenzofurans in a region in 
Spain containing an incinerator; Elliott et at^ studied 
cancer of the larynx and lung in Britain; Vinceti et aH 
focused on congenital abnormalities in Italy. None of 
these studies reported an effect. In a review of studies 
concerned with waste management, Rushton^^ states 
that the body of literature does not support concerns 
about proximity to waste disposal, since results vary for 
different cancer sites and many studies are unable to 
take into account confounding factors such as socio- 
economic status. 

The aim of this study was to investigate whether 
large-scale incinerators in England have an impact on 
the health of the people living close to them. Our study 
used a novel design in which circular case regions con- 
taining an incinerator were matched to similar control 
regions without an incinerator, rather than using the 
more common method of matching individuals. 
The main objectives were (1) to investigate whether 
there are differences in the incidences of health out- 
comes between case and control circles, taking into 
account some possible confounding factors; (2) to inves- 
tigate whether incidence within circles varies according 
to proximity to an incinerator and (3) to set a baseline 
for future work that allows for prolonged lead times in 
adult carcinogenesis. To undertake this, we conducted a 
retrospective analysis of the spatial distribution of 
selected adverse health outcomes in the vicinity of 
large-scale incinerators in England. We compared health 
outcome counts at a population level in case and 
control regions, adjusting for deprivation, age and sex, 
and using distance from the incinerator as a proxy 
measure of exposure. 

METHODS 

The study period was 1998-2008. We included all munici- 
pal waste incinerators in England with a capacity in 
excess of 150 000 tonnes per annum (tpa) of waste mater- 
ial that were operational on or before 31 December 1998 
and for at least 5 years subsequently. This time period was 
chosen as the EU Waste Incineration Directive^^ came 
into force in December 2003 for new plants and for exist- 
ing plants in December 2005, allowing time for 
retrofitting equipment to meet the directive. This identi- 
fied seven incinerators, but we excluded two incinerators 
sited in London as they were less than 20 km apart, 
making it impossible to identify suitable matching 
control regions, as described below. 

The unit of geography used in this study is the Lower 
Layer Super Output Area (LSOA) of the UK census. 
There are 34 378 LSOAs in England and Wales, with a 
minimum population of 1000 and a mean population of 
1500. Case circles around each incinerator were defined 
as the set of all LSOAs whose geometric centroid was 
contained within a circle of radius 10 km centred on the 
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incinerator. Guided by previous literature, we chose 
this radius to be of sufficient size to include the area in 
which any health effects are likely to be observable, 
without being so large that effects would be masked. 

Control circles were selected as follows. For each Local 
Authority (LA) containing an incinerator, we identified 
the most similar LA using criteria defined by the Office 
for National Statistics (ONS) and explained in detail on 
its website. This measure was created by combining 
census data variables representing key characteristics of 
the area, including: demographic structure, household 
composition; housing; socioeconomic character; employ- 
ment; and industry sector. The centroid of each selected 
LA was determined using Maplnfo, and a circle of 
radius 10 km centred on this point was drawn. The 
control circle was defined as the set of all LSOAs whose 
centroids fell within the selected circle, provided that no 
incinerator with a capacity in excess of 50 000 tpa lay 
within it or within 10 km of its boundary. When this 
occurred, the second most similar match was used, 
which was necessary only for the Stoke-on-Trent/ 
Sunderland pairing. The five case circles with their 
matched control circles are mapped in figure 1. We 
determined the geometric centroid of each LSOA in 
each case or control region, and used this to establish 
the distance of each LSOA from the incinerator. 

From the ONS we obtained, for each LSOA, annual 
counts of the following health outcomes, chosen as those 
most likely to have a biologically plausible association 
with incinerator by-products^^: childhood cancer inci- 
dence (<15 years, ICD-10 codes COO-97^^); childhood 



Case and Control Regions 




leukaemia incidence (<15 years, C91-95); leukaemia inci- 
dence (C91-95); liver cancer incidence (C22); lung 
cancer incidence (C34); non-Hodgkin's lymphoma inci- 
dence (C82-85); all-cause mortality; infant mortality 
(<1 year) and liver cancer mortality (C22). 

We also obtained LSOA-level population counts, by 
5-year age groups and sex, for 2001.^^ We used the 
English Indices of Deprivation 2004 as an index of 
multiple deprivation (IMD); this also uses 2001 census 
data. The IMD is a single summary measure based on 
seven domains of deprivation: Income Deprivation; 
Employment Deprivation; Health Deprivation and 
Disability; Education, Skills and Training Deprivation; 
Barriers to Housing and Services; Crime; and Living 
Environment Deprivation. 

Models 

We calculated expected counts for each LSOA in the 
study and for each health outcome using indirect stand- 
ardisation by applying the national age and sex-specific 
incidence or mortality rate to the age-sex structure of 
the LSOA.^^ We then considered two models. Model 1 
investigates the temporal trend at circle level and estab- 
lishes the importance or otherwise of including the tem- 
poral effect in subsequent analysis; Model 2 uses 
time-aggregated data and includes a function that allows 
for smoothly varying risk with increasing distance from 
the incinerator site. 

Model 1: temporal analysis 

Here, we use health outcome data aggregated over all 
LSOAs within each circle. The expected count obtained 
using indirect standardisation is included in the model 
as an offset, Eyt for the ith pair (i=l,...,5) and jth circle 
(j=l,2) in year t (t=l,...,ll), Wy represents IMD, Ci indi- 
cates membership of pair i and Vj indicates case/ control 
status (1 for case circles, 0 for control circles). The 
model assumes that the count, Yyt, for the ith pair and 
jth circle in year t is such that: 



Yijt Poisson(/XijJ 



where 



fi--^ = Eijtexp(Q;+^Wij + £i + i7Vj + Tt) 
and where - means 'is distributed as'. 




Case 


Control 


+ 1 Billingham 


A 1 Doncaster 


+ 2 Tyseley 


A 2 Wolverhampton 


+ 3 Coventry 


A 3 Derby 


+ 4 Eastcrofl 


A 4 Manchester 


+ 5 Stoke 


A 5 Sunderland 



Figure 1 Locations of case and control regions. 



IVIodel 2: multiplicative factor for distance 

We adapt the model used by Diggle et to explore the 
possible effect of proximity to the incinerator. As popula- 
tion and IMD data originate from the 2001 census, we 
are unable to include LSOA-level offsets corresponding 
to individual years in this model. We therefore aggregate 
over time: Yy^ denotes the total number of cases over 
the study period, in the ith pair, jth circle and kth 
LSOA, and Ey^ the corresponding expected count. 



Reeve NF, Fanshawe TR, Keegan TJ, etal. BMJ Open 2013;3:e001847. doi:1 0.1 136/bmjopen-201 2-001 847 



3 



Spatial analysis of health effects of large industrial incinerators in England 



Distance to the incinerator, dyk, is included as a covari- 
ate in a continuous function g(y,^,dijk) with two para- 
meters: one for the elevation in disease risk at the 
source (Q and one for the decay in risk as distance 
increases (y) . The model is 

Yijk Poisson(/Xijk), 

where 

Mijk = Eijk exp ( a + /3wijk + Si + T7Vj ) g( y, ^, dijk ) 

and 

g(y,^,dijk) = l + ^exp[-(dijk/y) ]vj 

We also fit this model with (dyk/y)^ replaced by dyk/y 
to check the effect of the functional form of distance on 
the results. 

Overdispersion 

Where overdispersion is observed (ie, variance greater 
than the mean), we use the method described by 
McCullagh and Nelder,^^ which accounts for this by 
adjusting the model-based SEs. As this method is only 
valid when there is no residual spatial correlation, we 
apply the diagnostic test for residual spatial correlation 
used by Fanshawe et al^^ comparing the sample correl- 
ation of the distances between LSOA centroids and the 
squared differences between their residuals to that 
obtained by repeated random assignment of residuals. 
We also carry out tests of deviance to compare different 
models, allowing for overdispersion.^^ 



Confidence intervals 

All CIs presented are 95% CIs, calculated as point esti- 
mates ±1.96 SEs. 

RESULTS 

In this section, we first describe the results of an initial, 
exploratory analysis and then report the results of the 
models outlined in the Methods section. 

Initial analysis 

A simple illustration of the potential effect of distance 
from the incinerator, and of any differences between 
case circles, is given by the cumulative standardised mor- 
tality/incidence ratio (SMR/SIR) for each health 
outcome for each circle as distance from the incinerator 
increases. This is shown in figure 2. This figure shows 
substantially higher SMRs/SIRs in some circles than in 
others, but does not suggest a trend in relation to dis- 
tance from the incinerator for most case circles for most 
health outcomes. SIR estimates for childhood leukaemia 
are subject to large sampling variability, especially at 
small distances, because counts are smaller than for the 
other outcomes. 

Table 1 summarises deprivation and the age-sex struc- 
ture for all case and control circles. Several variables 
used in the matching algorithm are also shown — for 
these variables, values are taken from the 2001 Census 
for either the LA in which the incinerator is located or 
the LA with which it is matched. 

Differences in IMD percentiles between matched 
circles are small. Across the whole study region, the 
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Figure 2 Estimated standardised mortality/incidence ratios (SIVIRs/SIRs) for each health outcome with respect to distance to the 
centre of the circle (case circles only). Line numbers correspond to case circles shown in figure 1. Vertical axis is plotted on 
logarithmic scale. 



4 



Reeve NF, Fanshawe TR, Keegan TJ, etai BMJ Open 2013;3:e001847. doi:1 0.11 36/bmjopen-201 2-001 847 



Spatial analysis of health effects of large industrial incinerators in England 



Table 1 Comparison of demographic factors in case and control circles 




Totals 






Pair 1 




Pair 2 




Pair 3 




Pair 4 




Pair 5 




m 


Case 


Control 


All 
























circles 


circles 


circles 


Case 


Control 


Case 


Control 


Case 


Control 


Case 


Control 


Case 


Control 


IMD (percentiles) 




























25th 


15.58 


15.86 


15.74 


16.13 


16.06 


18.92 


15.58 


13.14 


8.81 


12.26 


16.64 


16.07 


18.36 


50th 


28.79 


29.53 


29.24 


33.08 


28.34 


35.76 


31.18 


21.90 


20.33 


23.53 


29.98 


26.58 


30.94 


75th 


47.40 


45.62 


46.59 


53.20 


45.52 


51.13 


42.70 


35.71 


33.75 


42.88 


51.54 


42.72 


44.62 


Male (%) 


49.05 


48.96 


49.00 


48.52 


48.89 


48.93 


49.00 


49.76 


49.17 


49.43 


49.06 


48.63 


48.55 


Age 




























Mean 


41.40 


41.57 


41.48 


41.19 


41.46 


41.11 


41.53 


41.54 


41.65 


41.82 


41.48 


41.68 


41.82 


SD 


12.04 


12.05 


12.05 


12.00 


11.98 


11.95 


11.91 


12.14 


12.10 


12.30 


12.14 


11.93 


12.11 


Population (/100 000) 


26.57 


25.72 


52.29 


3.89 


2.25 


10.18 


5.89 


3.75 


2.83 


5.22 


10.48 


3.53 


4.27 


Population density 


27.45 


24.39 


25.92 


8.75 


5.05 


36.49 


34.07 


30.50 


28.41 


35.78 


33.97 


25.75 


20.43 


(people per hectare) 




























Not born in UK (%) 


9.23 


7.97 


8.60 


3.00 


2.87 


16.50 


11.24 


13.01 


8.55 


9.94 


14.78 


3.72 


2.40 


Single-person 


17.84 


16.80 


17.32 


14.56 


12.69 


18.65 


16.20 


16.52 


16.48 


23.39 


24.41 


16.09 


14.23 


household (%) 




























HE qualifications (%) 


15.20 


15.31 


15.25 


15.17 


11.78 


16.62 


13.63 


16.73 


17.73 


17.62 


21.42 


9.87 


11.97 


Unemployment (%) 


4.78 


4.67 


4.72 


4.98 


4.20 


5.71 


5.32 


3.96 


4.00 


5.25 


5.02 


3.99 


4.81 





mean age is 41.48 years, and 49% of all people in the 
study are male. There is little difference in the values of 
these factors averaged over all case circles, compared 
with all control circles. In each individual circle, the 
average age and the sex ratio are very similar to the 
overall average. Population counts between case and 
control circles within pairs vary, but we adjust for this 
using standardisation to calculate expected counts, as 
described in the Methods section. 

The first part of table 2 shows average unstandardised 
yearly numbers of cases for each outcome in both case 
and control circles. Most outcomes have similar average 
counts in case and control circles; apparent exceptions 
are leukaemia incidence and infant mortality, which 
have higher average counts in case circles (329 and 260, 
respectively) than in control circles (289 and 202, 
respectively), and lung cancer incidence, which is 
higher in control circles (2077) than in case circles 
(1861). Figure 3 shows temporal trends in SMRs/SIRs. 
Although all-cause mortality shows a gradual decline in 
raw count over time, its SMR shows an increasing tem- 
poral trend in both case and control regions. Similarly, 
increasing temporal trends in SMR/ SIR are apparent for 
infant mortality and lung cancer incidence. Temporal 
trends in case and control circles are similar for all 
outcomes. 



Model-specific analysis 

Model 1: temporal analysis 

The effect of year was statistically significant, showing an 
increasing trend over time, for leukaemia, liver cancer and 
lung cancer incidence, and all-cause, infant and liver 
cancer mortality (table 2, last three columns). For the out- 
comes showing a larger temporal effect, there were only 
small differences in the trend for case circles compared 
with control circles (figure 3). In addition, the population 
and IMD data originated from the 2001 census, and we 
were therefore unable to include the correct LSOA-level 
offsets corresponding to individual years in this model. For 



these reasons, we aggregated data over time before fitting 
models for spatial variation within circles. 

Model 2: multiplicative factor for distance 

The effect of distance from an incinerator for each 
outcome is indicated by the parameters for the elevation 
in risk at the source and the decay in risk (^ and y) . As the 
estimate of disease risk at a given distance is a combination 
of these two parameters, it is not possible to convert the 
individual parameter estimates to a simple relative risk 
(RR). Instead, the estimates of these two parameters are 
shown in table 3. Parameter estimates for non-Hodgkin's 
lymphoma incidence and liver cancer mortality had very 
large SEs and so are not reported. In addition, the likeli- 
hood function for liver cancer incidence was flat and 
showed no appreciable improvement in fit compared to a 
model with the g(y,^,dijk) term removed, so results from 
the latter model are shown for this variable (table 3). 
Table 4 shows the results of analysis of deviance for model 
comparison. There was no difference in fit between using 
squared distance and distance. Although there is some 
improvement in fit for infant mortality and lung cancer 
mortality for the model that includes distance compared 
to the model that ignores distance, there was no evidence 
of increased risk close to the incinerator for these or any 
other outcomes, as the ^ parameter that represents this 
increase in risk was either negative or not significantly dif- 
ferent from zero. 

There was a statistically significant increase in risk for each 
unit increase in IMD of 1.7% for infant mortality, 1.5% for 
lung cancer incidence and less than 1% for liver cancer inci- 
dence and mortality and all-cause mortality (table 3). 

The RR of living in a case circle compared with living 
in a control circle for leukaemia incidence was 1.137 
(95% CI 1.065 to 1.2 13), but it was not significant for any 
other outcome except all cause mortality, for which the 
effect is in a counter-intuitive direction, with RR 0.965 
(95% CI 0.947 to 0.984). For leukaemia incidence, there 
was no evidence of a gradient in risk according to dis- 
tance from the incinerator (table 3) . 
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Table 2 Total counts of health outcomes and relative risks (95% CIs) from Model 1 : Temporal analysis 



Health outcome 
9 


Total count 
All Case 
circles circles 


Control 
circles 


IMD (exp(p)) 




Case circle versus 
control circle (exp(Ti)) 


Year (exp(x)) 




v^riMunoou Cancer 


135 


70 


DO 




\ .\j \o) 


\ .UH-U ^u.yoD lO \ .\Ol ) 


n QQO (Ci on A \r\ 

u.yyu ^u.y/H- lo 


1 .UUDJ 


i n i H Q n r'Q 
11 lUlUUi lUc; 


















r^hilrlhr»nrl 
V^ilMUilUUU 


45 


23 


99 

C.C. 




1 ^9'^^ 

1 .UtlO) 


1 CiCiR (Ci P'^7 tn 1 9^R^ 
1 .UUO yJ.OOl lU 1 .^UDj 


1 nnzL ir\ Q7R tn 




IcUKaclllla 


















iiiciuci ice 


















1— UUlxctUl 1 Met 


618 


329 




n QQ7 (C\ QR7 tn 


1 ()n7\ 


1 . 1 OO ^ 1 .UOU lU I.I <J^) 


1 ni n M 009 tn 

1 .U 1 U \ \ X)\Jc. lU 


1 01 

1 .U 1 


inpirlonpo 

11 lUIUCI lUC 


















Liver cancer 


293 


142 


\ O 1 


\ .KjOc. ^ I .U I / lO 


1 .UH-Oj 


n QP9 (Ci Q1 Q tr» 1 nc;^^\ 


1 09A M 01 9 tr» 
1 .U^H- ^ 1 .U 1 ^ lO 


1 .UoOj 




















l-Uiiy Udiluei 


3939 


1861 


9077 


1 097 M 099 tn 


1 .UvJ 1 ) 


n Qn^^ /^n rr7 tn n 09*^^ 


1 01 1 M OOR tn 
1 .U 1 1 ^1 .UUO lU 


1 OIJ.^ 
1 .U 1 H-j 


iiiciuei ice 




















820 


406 


413 


w . ceo \\J . CUt^ IKJ 


1 .002) 




1 004 (Vi QQ7 tn 


1 .010) 


lymphoma 


















Incidence 


















All-cause mortality 


54009 


26750 


27259 


1.006 (1.005 to 


1 .007) 


0.958 (0.953 to 0.963) 


1.011 (1.010 to 


1.012) 


Infant mortality 


462 


260 


202 


1.015 (1.003 to 


1 .028) 


1.118 (1.054 to 1.186) 


1.025 (1.01 6 to 


1 .034) 


Liver cancer 


256 


129 


127 


1.035 (1.01 8 to 


1 .052) 


1 .008 (0.933 to 1 .089) 


1.014 (1.003 to 


1 .026) 


mortality 


















Estimates correspond to a unit change in the given variable. 



IMD, index of multiple deprivation. 



After adjusting for overdispersion, the fit was 
adequate. The residual spatial correlation was significant 
for childhood cancer and lung cancer incidence and 
infant mortality, with p values of 0.02, <0.001 and 0.01, 



respectively. For these outcomes, the correlations were 
all less than 0.019. As these correlations are small, we 
conclude that the spatial dependence is negligible, 
making the correction for overdispersion reliable. In 




Figure 3 Estinnated standardised mortality/incidence ratios (SMRs/SIRs) over time for each health outcome, averaged across 
case and control circles. Vertical axis is plotted on the logarithmic scale. 
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Table 3 Relative risks and parameter estimates (95% CIs) from IVIodel 2 



Health outcome IMD (exp(p)) 



Case circle versus 
control circle (exp(ri)) 



Change in risk at 
source (Q 



Rate of decay 
(km (7)) 



Childhood 0.999 (0.997 to 1 .002) 

cancer incidence 

Childhood 0.998 (0.993 to 1 .002) 

leukaemia 
incidence 
Leukaemia 
incidence 
Liver cancer 
incidence 
Lung cancer 
incidence 

Non-Hodgkin's 0.999 (0.998 to 1 .000) 
lymphoma 
incidence 
All-cause 
mortality 

Infant mortality 1 .01 7 (1 .01 5 to 1 .01 9) 
Liver cancer 1 .008 (1 .006 to 1 .01 0) 
mortality 



0.999 (0.998 to 1 .000) 
1.009 (1.007 to 1.010) 
1.015 (1.014 to 1.015) 



1 .059 (0.797 to 1 .407) 
1 .023 (0.704 to 1 .486) 

1.137 (1.065 to 1.213) 
0.984 (0.914 to 1.060) 
0.936 (0.874 to 1 .003) 
0.987 (0.945 to 1.031) 



1 .008 (1 .008 to 1 .009) 0.965 (0.947 to 0.984) 



0.974 (0.600 to 1 .582) 
1 .004 (0.928 to 1 .086) 



-0.281 (-4.404 to 3.842) 1 .715 (0 to 8.151) 

-0.650 (-12.320 to 11.020) 1.148 (0 to 7.120) 

-0.472 (-6.043 to 5.099) 0.682 (0 to 2.891 ) 

-0.096 (-0.1 88 to -0.004) 3.598 (0 to 8.624) 



0.187 (-0.469 to 0.843) 
0.225 (-0.529 to 0.980) 



0.762 (0 to 1.944) 
6.454 (0 to 15.715) 



For outcomes where no estimates are presented, see text for details. Relative risks correspond to a unit change in the given variable. 



addition, plots of residuals (pooled across all circles) 
against distance from the incinerator and against direc- 
tion from incinerator did not reveal any systematic 
trends, suggesting that the modelling assumptions are 
reasonable (see online supplementary material). 

DISCUSSION 

We have presented the results of a study of the health 
effects of large industrial incinerators in England 
between 1998 and 2008. Across a broad range of out- 
comes that have been linked to pollutants released by 
incinerators in previous studies, our results show no ele- 
vated risk for individuals living in areas containing an 
incinerator compared to individuals living in matched 
areas without an incinerator. 

In the models that constitute the main analysis, we 
made allowances for both the age-sex structure of the 
underlying population, through indirect standardisation, 
and the key confounder of deprivation, measured at the 
small-area (LSOA) level. The results of fitting a continu- 
ous risk function to reflect the patterns of risk of specific 
outcomes across regions containing an incinerator need 
to be interpreted carefully. There was an apparently 
increased risk of leukaemia incidence for those living in 
case regions, but for this outcome the distance effect was 
non-significant and the trend, counter-intuitively, was 
towards a decreased risk in the vicinity of the incinerator. 

For all outcomes considered, within case circles, there 
is little evidence of an increase in risk for those living in 
proximity to an incinerator compared to those living 
further away. One might have anticipated a dose- 
response relationship with distance from the incinerator 



if the incinerators were responsible for causing an eleva- 
tion in disease or mortality risk. The fact that no such 
relationships emerged supports our general conclusions. 

Others have used dispersion modelling of 

emissions to further define exposure, particularly of 
dioxins, proven carcinogens. In this study, we chose to 
use a distance-based proxy measure of exposure because 
we did not have accurate retrospective data on location- 
specific emissions, and there remain doubts about the 
accuracy of air dispersion model outputs in this setting, 
especially when they cannot be validated against reliable 
measurements from monitoring networks. 

Part of the background to this study is the continuing 
community concern about incinerators, a concern 
which does not discriminate incinerators by emissions. 
Although it is possible that incinerator emissions differ 
between sites, they are likely to be more consistent over 
time and space than is sometimes suggested, combustion 
conditions being more important than precursor waste 
content. This supports our decision to treat incinerators 
at different sites as if they are homogeneous units, 
although differences in emissions between sites or over 
time may reduce the power of our study to identify 
authentic health effects. 

Although some of the cancers we examined may have 
longer lead times from exposure to clinical disease than 
the period covered by the study, infant mortality and 
some of the childhood cancers have much shorter lead 
times, less than the 10 years covered by the period of 
the present study. For the adult cancers, repeating our 
study at an appropriate interval, based on the same sites 
which are now subject to more stringent controls than 
before 2005, may go some way to address this issue. 
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Exploratory analysis showed substantial differences in 
case numbers between different circles, possibly due to 
factors such as increased levels of deprivation. This 
emphasised the importance of including the pairing as a 
factor in the model. While no clear effect of distance 
from the incinerator was observed, it is possible that a 
distance effect was masked by other factors not taken 
into account in these plots, such as the deprivation of 
the area. 

Higher counts in case circles than control circles were 
observed for leukaemia incidence and infant mortality 
(table 1). This may indicate a potential incinerator 
effect, but may also be at least partially explained by 
other factors. In particular, this initial analysis makes no 
adjustment for the age-sex structure of the regions or 
differences in deprivation level. 

Several other factors relating to the design and nature 
of the study may explain our findings. In common with 
other studies in this field,^ ours may suffer from unob- 
served confounding factors that our analysis has been 
unable to address. Although the dataset, based on offi- 
cial data from the ONS, is as comprehensive as could be 
expected for a study of this type, our analysis is based on 
small-area-level data, and therefore we are unable to 
adjust for individual-level risk factors such as smoking 
status, diet and genetic risk factors, which are strongly 
associated with many of the outcomes considered.^^"^^ 

Higher levels of deprivation are strikingly linked to 
increased incidence in many of the disease outcomes we 
have considered: all-cause and infant mortality, lung 
cancer incidence and liver cancer incidence and mortal- 
ity. The decision to employ a circle-matched design in 
this study, in which case and control circles were also 
matched on a variety of other socioeconomic factors, 
reduced the need to adjust for many of these factors 
explicitly in our analysis but does not entirely eliminate 
it. Also, while the IMD is a generally accepted summary 
measure of deprivation, it cannot fully capture aspects 
of deprivation that may affect the risk for specific dis- 
eases. In addition, we did not take into account other 
sources of local air pollution. This may be a beneficial 
avenue for future studies to explore. 

The statistical methods used in our analysis include a 
continuously decaying distance function model, which 
provided an adequate fit and gives additional informa- 
tion about the likely effect of distance on disease risk. 
Diggle et and Lawson^^ have suggested fitting step 
function models, which include a plateau of constant 
elevated risk up to an unknown distance 5, and models 
allowing for angular direction from a point source. 
When we attempted to fit these models to our data, we 
found that parameter estimates were unstable. Further 
research is needed to evaluate the practical feasibility of 
fitting these models to area-level data. 

In summary, our findings broadly support those of 
previous research,^ even though the residential his- 
tories of the individuals contributing to this study are 
unknown and many of the cancers considered have long 
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lead times. It may require a carefully constructed 
individual-level cohort study, possibly also using validated 
dispersion modelling of emissions, to address these con- 
cerns definitively 
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